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The construction of initial-data sets representing binary black-hole configurations in quasi-circular orbits is studied in the 
context of the conformal-imaging formalism. An effective-potential approach for locating quasi-circular orbits is outlined for 
the general case of two holes of arbitrary size and with arbitrary spins. Such orbits are explicitly determined for the case of two 
equal-sized nonrotating holes, and the innermost stable quasi-circular orbit is located. The characteristics of this innermost 
orbit are compared to previous estimates for it, and the entire sequence of quasi-circular orbits is compared to results from the 
post-Newtonian approximation. Some aspects of the numerical evolution of such data sets are explored. 
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I. INTRODUCTION 

The numerical simulation of binary black- hole systems 
begins with the specification of appropriate initial data. 
A general method for specifying the initial data of bi- 
nary black-hole systems has been described in Cook et 
al. Q (hereafter Paper I). Because of the circularizing 
effects of gravitational radiation damping, we expect the 
orbits of most tight binary systems to have small eccen- 
tricities. Therefore, a method is needed which can dis- 
cern which data sets, within the very large parameter 
space of binary black-hole initial-data sets, correspond 
to binary black holes in a quasi-circular orbit. In this pa- 
per, I will describe such a method and explicitly compute 
the initial-data parameters necessary for describing the 
quasi-circular orbit of two equal-sized nonrotating black 
holes. In addition, this method yields an estimate of the 
innermost stable quasi-circular orbit for two equal-sized 
black holes. 

The general framework being used for defining initial- 
data sets containing black holes is known as the 
conformal-imaging formalism 0M0,E1@M|. It is based 
on the Arnowitt-Deser-Misner (ADM) (8|, or 3+1, de- 
composition of Einstein's equations, York's conformal 
and transverse-traceless decompositions of the constraint 
equations, and a method of imaging applicable to tensors. 
Application of this approach to the case of two black holes 
with arbitrary linear and angular momenta on each hole 
has been explored in Paper I. In that work, three indepen- 
dent numerical approaches were described for construct- 
ing initial-data sets. I will make use, in this paper, of 
the "Cadez-coordinate approach" for solving the Hamil- 
tonian constraint, a three-dimensional (3D) quasi-linear 
elliptic partial differential equation. 

Using numerical initial-data sets, an "effective poten- 
tial" can be constructed which consists of the gravita- 
tional binding energy between the holes plus the kinetic 
energy of the holes, with the restriction that all physical 
parameters characterizing the system are held fixed ex- 



cept for the separation of the holes. If we fix the momenta 
of the holes appropriately, then an initial-data set occur- 
ring at a minimum of the effective potential will represent 
two black holes in a quasi-circular orbit. For the case 
of two equal-sized holes with no intrinsic rotation, the 
set of physical parameters that must be held fixed con- 
tains only the orbital angular momentum of the system. 
For this simple case, we have located those initial-data 
sets which represent two black holes in quasi-circular or- 
bit at various separations. As expected P|l(|, there is 
some minimum separation required before a minimum of 
the effective potential exists. This minimum separation 
serves as an estimate for the point at which the secular 
inspiral of the two holes ends and the final dynamical 
plunge to coalescence occurs. 

I will begin with a description of the parameter space of 
3D initial-data sets, a brief description of how initial-data 
sets are constructed, and how the physical parameters as- 
sociated with a given data set are computed. I will con- 
tinue with a description of the effective-potential method 
as applied to initial data constructed via the conformal- 
imaging approach. Using this method, I explicitly locate 
the sequence of quasi-circular orbits (as a function of sep- 
aration) for equal-sized holes and compare this sequence 
to the results obtained via a post-Newtonian analysis. 
Also, the characteristics of the innermost stable quasi- 
circular orbit are compared to the results of previous es- 
timates for this orbit. I conclude with some observations 
about the eventual evolution of initial-data sets repre- 
senting two equal-sized holes in quasi-circular orbit. 



II. INITIAL DATA 

A detailed description of how binary black- hole initial- 
data sets are constructed within the conformal-imaging 
formalism is given in Paper I and references therein. In 
brief, the configuration is parameterized first by fixing 
the locations, Ci and C2, of the two holes in a Cartesian 
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flat conformal background space along with the radii, a± 
and «2) °f the two holes. If we let a\ set the fundamental 
length scale of the problem, then we can parameterize 
the physical locations and relative sizes of the holes in 
the background space by two dimensionless parameters, 
a and f3, defined by 



a 



ai 

a-2 

ICj-Cal 

ai 



(1) 
(2) 



In addition to the locations and sizes of the holes, we 
also specify for each hole both a linear momentum vector 
Pi ,2 and an angular momentum vector Si, 2. The physical 
meaning of these momenta is seen in the following way. 
For a single isolated hole, the physical linear and angular 
momenta contained in the initial-data set, as measured 
at infinity, are given directly by the parameters P and 
S. For two holes, the total physical linear momentum of 
the system is simply the vector sum Pi + P2. The total 
physical angular momentum is given by the vector sum 
Si + S2 + J, where J is the orbital angular momentum of 
the system which will be described in detail below. With 
a\ setting the fundamental length scale, we find that the 
momenta of two black holes within the initial-data set 
is completely specified by the four dimensionless vector 
parameters: Pi/ai, P2/0.1, Si/af, and S2/a 2 . 

Aside from the complications of a given numerical tech- 
nique used to solve for the initial data, we see that the 
specification of a general two-hole configuration requires 
that fourteen dimensionless parameters be fixed. Hav- 
ing chosen these fourteen parameters, an initial-data set 
is constructed by solving the momentum and Hamilto- 
nian constraints of general relativity as outlined in Pa- 
per I. The particular numerical technique used to solve 
the Hamiltonian constraint in this work is the "Cadez- 
coordinate approach" described in detail in Sec. IIIA of 
that paper and references therein. 

Once the initial-data set has been computed, various 
physical parameters characterizing the data set can be 
computed. In particular, the following quantities are cal- 
culated: the ADM energy of the system E, the proper 
surface areas of the marginally outer-trapped surfaces 
defining each individual hole Ai t 2, the dipole moment of 
the energy distribution d, and the shortest proper sepa- 
ration between the two marginally outer-trapped surfaces 
£. Definitions for the ADM energy and the dipole mo- 
ment can be found in Eqns. (24) and (25) of Ref. 0. 
Note that it is actually the dimensionless ratios E/a%, 
Ai ,2/0-1, d/af, and l/ai which are computed. These 
quantities, together with the initial-data parameters, al- 
low us to compute the orbital angular momentum J of the 
system. Generalizing the calculation found in Ref. O] , it 
is straightforward to show that the orbital angular mo- 
mentum for a configuration of two black holes is given 
by 



where O is the point in the background space about 
which the angular momentum is defined. The only 
unique choice for O is the center of energy of the sys- 
tem which, if chosen as the origin of coordinates, results 
in a vanishing dipole moment. Using the definition of the 
ADM energy and the dipole moment, it follows that the 
center of energy is at O = d/E. Therefore, the general 
definition of the orbital angular momentum is 



Ci-- ) xPi 



E 



x P, 



(4) 



III. THE EFFECTIVE-POTENTIAL METHOD 

The definitions of the ADM energy, total linear and 
angular momenta at infinity, the dipole moment, the 
proper separation of the holes, and the areas of the two 
marginally outer-trapped surfaces are rigorously defined 
physical quantities. To define the effective potential of a 
configuration, it is necessary to use the concepts of the 
masses and spins of the individual black holes and to have 
a measure of the effective binding energy between the two 
holes. However, none of these quantities are rigorously 
defined in general relativity for a strong-field nonstation- 
ary configuration. 

In the limit of large separations and small linear mo- 
menta and spin on the holes, the following definitions 
hold. The mass of each hole can be defined via the 



Christodoulou formula 12 



M 2 = Ml 



S 2 



4M ; 2 



(5) 



with the irreducible mass M iT « ^A/lGir and S be- 
ing the magnitude of the spin on the hole. As shown 
in the Appendix, we can approximate S for each hole 
by the magnitude of its respective spin parameter | S1.2 1 - 
Though it seems possible that the linear momentum on 
one hole could induce a spin on the other hole, this is, in 
fact, not the case. Thus, the quantity J, which we identi- 
fied above as the orbital angular momentum, is not con- 
taminated by an induced spin on the holes. We are there- 
fore justified in defining the orbital angular momentum 
of the system by equation (^), and we define the spins 
of the individual holes via their respective spin param- 
eters. Finally, the effective binding energy E^ between 
the holes is defined as 



E h = E - Mi - M 2 , 



(6) 



J = (Ci - O) x P x + (C 2 - O) x P 2) 



(3) 



where Mi ,2 are the masses of the two holes as defined 
above. Note that the effective binding energy contains 
both the gravitational binding energy between the two 
holes and their orbital kinetic energies, but not the rota- 
tional kinetic energy of the individual holes. 

These definitions are rigorous only in the limit of 
infinite separation and zero momenta on either hole. 
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The limit of zero momenta on the holes is required be- 
cause initial-data sets containing a single black hole con- 
structed via the conformal-imaging approach necessar- 
ily contain some spurious gravitational- wave energy [fl3| . 
The same is true of multihole initial-data sets, however, 
Cook and Abrahams Q] have shown that the spurious 
radiation content for the case of two holes is quite small 
so long as the holes are modestly separated and the mo- 
menta are not excessively large. We will find that these 
constraints are satisfied for the majority of configurations 
of physical interest. 

Henceforth, we will take the masses of the holes, their 
spins, mutual binding energy, and orbital angular mo- 
mentum to be defined as given above. 

We turn now to the definition of an effective potential 
useful for determining the location of quasi-circular or- 
bits. In general, such an effective potential should be a 
function of the separation and sizes of the holes, the or- 
bital angular momentum of the system, the spins of the 
holes, and the gravitational radiation content of the sys- 
tem. Within the conformal-imaging approach, one has no 
freedom in specifying the radiation content of the system. 
This is fixed by the demands that the spatial 3-metric be 
conformally flat and that all fields satisfy an isometry 
condition (cf. Ref . ) . With this restriction, we see that 
the effective potential should be a function of nine phys- 
ical parameters but a general two-hole initial-data set 
depends on fourteen initial-data parameters. 

To reduce the size of this parameter space, we first de- 
mand that the configuration be in a center of momentum 
frame. This restriction requires 



Pi 



0. 



(7) 



A configuration representing a quasi-circular orbit should 
satisfy 



Pi, 2 • (C2 — Ci) 



0. 



(8) 



Together, equations (Q) and (||) reduce the fourteen- 
dimensional initial-data parameter space to nine dimen- 
sions. These parameters are a, /3, the magnitude of the 
linear momentum on either hole P/ai, Si/af , and S2/CI?. 

The respective dimensionless physical parameters of 
the effective potential are X = M1/M2, i/m, J/fxm, 
Si /Mi, and S2 /M| where are the masses of the 

individual holes, m = M± + M2 is the total mass, 
fx = M\Mijm is the reduced mass, and J is the mag- 
nitude of the orbital angular momentum of the system. 
Finally, the dimensionless effective potential is given by 
the binding energy as E^/fj,. 

Finding quasi-circular orbits is now a conceptually easy 
task. We compute as a function of £/m while hold- 

ing the remaining eight physical parameters constant. A 
minimum in Eb/(J> then corresponds to a "stable" quasi- 
circular orbit. In addition to locating the quasi-circular 
orbits, we can also estimate the orbital angular velocity 
17 of the system as measured at infinity. This is given by 
taking the derivative of the binding energy with respect 



to the orbital angular momentum while holding all other 
parameters fixed. In dimensionless form then one obtains 



m 17 



dJ '/ ' \xm 



(9) 



Though conceptually straightforward, the computa- 
tional task of locating quasi-circular orbits is difficult. 
The main difficulty arises from the fact that the physical 
parameters that must be held fixed are not independently 
correlated with their respective initial-data parameters. 
That is, holding eight of the nine initial-data parame- 
ters fixed while varying f3 will not result in an effective- 
potential curve. We see then that, in general, the prob- 
lem of determining one quasi-circular orbit is quite in- 
volved. It requires finding the roots of eight functions, 
each of which depends on eight parameters, at each value 
of the separation at which the effective potential is eval- 
uated. 

Fortunately, the size of this parameter space can be 
cut in half. The only definition we have available for the 
direction of the spins of the holes (see the Appendix) is 
the direction of their respective initial-data spin param- 
eters. The directions of the spins are then fixed relative 
to the separation of the holes and the plane of the orbit, 
which arc defined by Pi, 2 and C2 — Ci. This direction is 
independent of the numerical solution of the Hamiltonian 
constraint, so only the magnitude of the two initial-data 
spin-vector parameters needs to be varied to hold the 
physical parameters fixed. As a result, for a fixed value 
of (3 the following four equations must be satisfied: 



X, 



J 



Si 



Ml 
S 2 



Ml 



(a.P/a 1 ,S 1 /a 2 1 ,S 2 /al;P) 
(a,P/ai,Si/a*,S 3 /al;P) 
(a,P/ai,Si/a*,Sa/al;P) 



X 
J 

Ml 
Ml 



(10a) 
(10b) 

(10c) 

(lOd) 



When these four equations are satisfied, equation (g) 
yields a value of the effective potential E^/n at some 
value of the physical separ ation l /m. Changing the value 
of j3 and resolving Eqns. (10a)-(10d) results in another 
value of the effective potential at a different separation. 



IV. EQUAL-MASS NONROTATING HOLES 

The simplest application of the effective-potential 
method is to the case of two equal-sized black holes with 
no intrinsic spin. In this problem, S\/a\ = Sija\ = 
and, because of the symmetry between the holes, we 
know that a = 1 — > X = 1. Therefore, solvi ng fo r the 
effective potential requires solving only Eqn. (10b) as a 
function of P/a\ alone for a given [3. 
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The method used to solve equation ( |lOq ) is the follow- 
ing. For a given value of /3, the initial- value equations are 
solved at a sufficiently large number of values of P/ a\ in 
order to encompass all of the values of J/pm at which we 
want to evaluate the effective potential. Using interpola- 
tion, we can estimate new values of P/a\ that will yield 
solutions near the desired values of J//im; the procedure 
is repeated until any errors introduced by interpolation 
are sufficiently small. The most difficult part of the prob- 
lem is to solve the initial-value equations with sufficient 
accuracy. Typically, both the ADM mass E/ai and the 
areas of the marginally outer-trapped surfaces Ai_ 2 /a 2 
need to be determined to a relative error of ~ 1CT 5 . 

Currently, the only numerical method capable of solv- 
ing the initial-value equations to this accuracy is the 
multigrid-based "Cadez-coordinate approach" described 
in detail in section IIIA of Paper I. Such high accuracy 
can be obtained through the use of Richardson extrapo- 
lation. As described in Paper I, the Cadez-coordinate ap- 
proach used to solve the Hamiltonian constraint results 
in a numerical solution for the conformal factor ip num , 
which has an asymptotic (h — > 0) expansion given by 



V> num = V + h 2 et + h 4 ef 



(11) 



where ip is the analytic solution of the Hamiltonian 
constraint, h is the basic scale of discretization, and 
are /i-independent functions. In addition, the 
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numerical integrals for E/ai, d/a 2 , Ai^/af, and £/oi 
have all been constructed to yield analogous error ex- 
pansions that depend strictly on powers of h 2 . 

One final source of error which must be examined 
comes from the necessity of imposing an approximate 
outer-boundary condition (cf. Ref. Q). In order to min- 
imize the effects of this approximation, the outer bound- 
ary has been placed at a distance of at least 2000ai from 
the holes. 

Figure |l| displays a set of effective-potential curves for 
a wide range of values for J / ' pm. The displayed curves 
are interpolated results derived from 3000 Richardson- 
extrapolated data points, each resulting from the ex- 
trapolation of three separate solutions of the initial- value 
equations at resolutions similar to those described in 
Paper I. All solutions were generated on an IBM-SP1 
parallel computer and required a total computational 
time in excess of 3000 CPU hours. The value of J/ pm 
is held fixed along each of the thin curves which plots 
the effective potential E-^/p as a function of the proper 
separation of the holes £/m. The bold curve crossing 
several of the effective-potential curves represents a se- 
quence of quasi-circular orbits. This can be seen more 
clearly in Fig. |^ where the region containing minima in 
the effective-potential curves is shown in expanded form. 
The bold line representing the sequence of quasi-circular 
orbits begins at the right at an l/m ~ 14. This line 
should, of course, extend to larger values of l/m, but 
data has not been computed in this regime. Following 
the sequence of quasi-circular orbits to smaller values of 



FIG. 1. The effective potential E^/p as a function of sep- 
aration l/m for the following values of the orbital angular 
momentum J/pm: 1.5, 2, 2.5, 2.75, 2.95, 2.976, 2.985, 3, 
3.05, 3.15, 3.25, 3.37, 3.5, 3.62, 3.75, 3.85, 4, 4.25, and 4.5. 
These values of J/pm label, respectively, curves from the bot- 
tom of the figure to the top. Also, plotted as a bold line is 
the sequence of quasi-circular orbits which cross the effec- 
tive-potential curves at local minima. 



FIG. 2. An enlargement of the section of Fig. [I] which con- 
tains the sequence of quasi-circular orbits. This sequence be- 
gins at the innermost stable quasi-circular orbit near l/m = 5 
on the J/pm = 2.976 curve and extends in the direction of 
larger separation. In the figure, the sequence terminates at 
the minimum of the J I pan = 3.85 effective-potential curve, 
although it should continue. 
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TABLE I. Physical and initial-data parameters charac- 
terizing certain configurations along the sequence of stable 
quasi-circular orbits. The data sets represented in this table 
have been constructed using the minus isometry condition of 
the conformal-imaging approach. 



e/m 




J/fim 




P/ai 




4.880 


-0.09030 


2.976 


0.172 


1.685 


11.82 


5.365 


-0.08890 


2.985 


0.145 


1.392 


13.28 


5.735 


-0.08684 


3.000 


0.130 


1.230 


14.43 


6.535 


-0.08112 


3.050 


0.104 


0.9868 


16.99 


7.700 


-0.07226 


3.150 


0.0774 


0.7752 


20.84 


8.695 


-0.06534 


3.250 


0.0622 


0.6613 


24.21 


9.800 


-0.05862 


3.370 


0.0504 


0.5734 


28.04 


10.96 


-0.05270 


3.500 


0.0414 


0.5066 


32.15 


12.02 


-0.04810 


3.620 


0.0352 


0.4609 


35.93 


13.16 


-0.04388 


3.750 


0.0300 


0.4218 


40.06 


14.07 


-0.04104 


3.850 


0.0270 


0.3960 


43.38 



i/m, we find that the minimum in the effective potential 
vanishes. At this point, the sequence of quasi-circular 
orbits terminates at an innermost stable orbit. 

For i/m < 4 we notice that some of the effective- 
potential curves pass through a local maximum and a 
second minimum. These additional minima should not 
be interpreted as a new sequence of stable quasi-circular 
orbits. Rather, this behavior indicates that the approxi- 
mations outlined in Sec. Ill, especially the identification 
of equation (|^) as a measure of the binding energy be- 
tween the holes, are breaking down. This assertion is jus- 
tified on the following grounds. Anninos et al. have 
shown that for time-symmetric initial data constructed 
via the conformal-imaging approach an event horizon en- 
compasses both holes on the initial-data slice when their 
separation parameter fx (not to be confused with the re- 
duced mass) is less than about 1.8. This corresponds to 
a separation parameter of (3 ~ 6.2 and to a proper sep- 
aration of t/m ~ 3. For the case of orbiting holes, we 
expect the effects of angular momentum to delay the on- 
set of formation of a common event horizon based on both 
physical intuition and the axisymmetric results of Cook 
and Abrahams H|. Thus, we should expect a common 
event horizon for t/m < 3 and, further, that the effective- 
potential method being used should not be trusted for 
t/m < 4. 

The characteristics of initial-data sets representing 
quasi-circular orbits of two equal-sized, nonrotating black 
holes are given in Table |[ In addition to the physical pa- 
rameters characterizing the system (t/m, E^/n, J/fim, 
and mSY), the table contains the initial-data parameters 
P/a± and j3 required to reproduce these particular data 
sets. Note that all initial-data sets described in this paper 
satisfy the minus isometry condition of the conformal- 
imaging approach. Values in this table should be con- 
sidered accurate to better than 1% with the exception of 
mf2, which should be considered accurate to a few per- 
cent. 



In order to gauge the accuracy with which we have lo- 
cated the innermost stable quasi-circular orbit, note that 
in the limit of a test mass orbiting a Schwarzschild black 
hole, the proper separation between the event horizon 
and the test mass is found to be 



— = 2 In 

m 



4.58. 



(12) 



This should be compared with a value of i/m — 4.88 
obtained from Table ffl. The ratio of these two values 
is 0.94. Kidder et al. |l(J obtain an analogous ratio of 
0.96, where separation is measured in terms of harmonic 
or deDonder coordinates. Kidder et al. obtain this re- 
sult via a critical point analysis of the equations of mo- 
tion through (post) 2 -Newtonian order. After altering the 
equations of motion significantly to reproduce exactly the 
test mass limit, Kidder et al. find that the ratio has 
changed to 0.83 and that the innermost stable circular 
orbit is characterized by the following physical parame- 
ters: E h /n 0.0378, J / urn ~ 3.83, and mfl ~ 0.0605. 

A comparison of these values with Table | shows that 
Kidder et al. are finding an innermost stable circular or- 
bit in which the holes are much farther apart than we 
find with the effective-potential method. In contrast to 
this, Blackburn and Detweiler |J have used a variational 
principle together with the assumption of a periodic so- 
lution to Einstein's equations to obtain an estimate for 
the innermost orbit for two equal-sized holes. Using a 
single trial geometry, which they call "rather unsophisti- 
cated" , they obtain an innermost orbit characterized by 
Eh/fJ' ~ —0.65, J//im ~ 0.85, and mfl ~ 2. Such an or- 
bit is much more tightly bound than seems possible from 
the effective-potential method. Blackburn and Detweiler 
point out that the assumptions of their variational prin- 
ciple have been violated by the time this innermost circu- 
lar orbit is reached. However, they describe a less tightly 
bound circular orbit that should not be in violation of 
the underlying assumptions of their approach. This or- 
bit has a binding energy of E\,//j, ~ —0.28, which is still 
three times larger than the binding energy obtained in 
this paper for the innermost stable quasi-circular orbit. 



V. THE POST-NEWTONIAN LIMIT 

Comparison with previous estimates for the innermost 
stable quasi-circular orbit of two equal-mass nonrotat- 
ing black holes is far from yielding a consensus as to its 
proper value. However, we can gain some insight into the 
reliability of the results derived in this paper by com- 
paring the sequence of quasi-circular orbits against the 
post-Newtonian description of circular orbits. Based on 
the (post) 2 -Newtonian results of Kidder et al. [jl6| for 
the binding energy, angular momentum, and equations 
of motion for a binary system with a circular orbit, it 
is straightforward to show that the binding energy and 
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FIG. 3. The effective potential Eb/fi as a funcion of the or- 
bital angular momentum J/fim for quasi-circular orbits. The 
solid line corresponds to the sequence of quasi-circular orbits 
computed in this paper. The long dashed line is the result 
obtained from Newtonian theory. The short dashed line is 
the result based on (post) ^Newtonian theory, and the dotted 
line is the result based on (post) 2 -Newtonian theory. 



angular momentum of two compact objects in a circular 
orbit must satisfy 
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(mQ) 



-2/3 



l + -(9 + rj) (mfi) 



2/3 



\4/3 



(15) 



where 77 = [i/m. The three terms inside the square brack- 
ets represent the Newtonian results for circular orbits 
along with the first and second post-Newtonian correc- 
tions. 

Figure [| compares the numerical results for binding 
energy versus orbital angular momentum for the se- 
quence of quasi-circular orbits against equation ( [13] ) with 
rj = 1/4. The numerical data is displayed as a bold solid 
line with cross marks denoting actual data points. The 



FIG. 4. The orbital angular momentum J//xm as a function 
of the orbital angular velocity mQ. Lines are as indicated in 
Fig. |. 



long dashed line corresponds to the Newtonian result, the 
short dashed line to the Newtonian result together with 
the (post) 1 -Newtonian corrections, and the dotted line to 
the full (post) 2 -Newtonian result. Notice that for large 
J I [im (large separation), the post-Newtonian expansion 
appears to be converging quite well toward the numerical 
result. Figures ^ and]^ are analogous plots for the orbital 
angular momentum versus orbital angular frequency and 
binding energy versus orbital angular frequency, respec- 
tively. Again, in the limit of large separation (small 
mf2) , the post-Newtonian expansions appear to converge 
well with the numerical results. This agreement at large 
separations, together with the general agreement in the 
shapes of the curves when the separation decreases, con- 
firm that the basic premises adopted in the definition of 
the effective-potential method are sound. For small sep- 
arations, we know that the approximations n eede d to de- 
fine the effective potential as outlined in Sec. Ill must be 
questioned. Unfortunately, the post-Newtonian approx- 
imation also breaks down as the separation of the holes 
becomes small and, it is impossible to gauge the valid- 
ity of the approximations based on the post-Newtonian 
results. 



VI. DISCUSSION 

The most important use of the results obtained in this 
paper is to narrow the range of initial-data parameters 
which must be considered in setting up an actual numer- 
ical simulation of the inspiral and collision of two black 
holes. Assuming one is interested in evolving initial data 
that represents something similar to a quasi-circular or- 
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FIG. 5. The effective potential Eh/fi as a function of the 
orbital angular velocity mD.. Lines are as indicated in Fig. 



the resulting black hole and that half of the angular mo- 
mentum is, nevertheless, radiated away during the final 
plunge and coalescence, we find that Jf/M f 2 ~ 0.4. We 
see then that a black hole resulting from the inspiral and 
final plunge of two nonrotating black holes must certainly 
be considered to be rapidly rotating, but it will not vio- 
late the Kerr limit. 

Now consider the minimal requirements of a numerical 
evolution of binary coalescence in the case of two equal- 
sized, nonrotating black holes. If the simulation is to 
remotely resemble the final plunge of two holes following 
a secular inspiral, then our best guess at initial data is 
that for the innermost stable quasi-circular orbit. From 
the estimate of the orbital angular velocity in Table |, we 
find that the orbital period is r ~ 37m. It is reasonable 
to assume that the final plunge will occur on a time scale 
comparable to that of the innermost orbit. If this initial 
data leads to an evolution that begins at the most one 
orbit before the beginning of the final plunge and we add 
to the evolution time a period sufficient to watch some of 
the ring-down, then we find that the numerical simulation 
must be capable of evolving for 90-130m. 



bit, then the results presented in Table Q allow us to esti- 
mate the minimum length of time that the system must 
be evolved, as well as some limits on the final state of the 
system. 

Beginning with the latter, let us assume that the bi- 
nary system coalesces and settles down to a Kerr hole 
with mass Mf and angular momentum Jf. From equa- 
tion (||) we know that the Kerr ratio Jf /Mf is 



A 

Mf 



(16) 



where we now let M- 1T denote the irreducible mass of 
the final Kerr hole. Assuming only that the individual 
holes do not rotate, we know that the irreducible mass is 
bounded by 



Mf + Ml 



(17) 



where Mi, Mi, etc. arc defined as before. Since we know 
that Jf < J, we find that 



Jf < V 

Ml ~ 1 



J 



2-q [im 



(18) 
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APPENDIX: THE MOMENTA OF INDIVIDUAL 
HOLES 

A rigorous quasi-local definition of momentum requires 
the presence of a Killing vector field £i s . The magnitude 
of the momentum associated with £^ within a given sur- 
face contained in a spatial slice is denoted ELm , where 



II(fc) = HiC(k) - 



Kl-8lK)? (k) d 2 S 3 . (Al) 



Examining equations (|l6|) and ( |l8|) we find that, so long 
as Jf/Mf < 2, the Kerr ratio Jf/Mf is guaranteed to Here, K\ is the extrinsic curvature of the slice and K 



be less than unity. For equal-sized holes, this implies 
that the Kerr ratio is satisfied as long as J / [im < 4. 
If we evolve initial data representing the innermost sta- 
ble quasi-circular orbit, we find that Jf/Mf < 0.958. If 
we make the severe assumptions that all of the bind- 
ing energy "released" in the coalescence is recaptured by 



is its trace. A general spatial slice contains no Killing 
vectors and only the total momenta of an asymptotically 
flat spatial slice can be defined. This definition requires 
integrating equation (|Al| ) at spatial infinity using asymp- 
totic Killing vectors £5 fc ^ , which are Killing vectors of the 
flat metric to which the physical metric is asymptotic. 
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Note, however, that the angular momentum will only be 
gauge-invariant if is an exact symmetry of the phys- 
ical metric (cf York jlTj ) . 

Now consider the momenta of gravitational initial data 
constructed within the conformal-imaging approach. Fol- 
lowing this approach, we conformally decompose the 
physical metric of a spatial slice as 7^ = ip^Jij ) where jij 
is the flat conformal background metric and tp is the con- 
formal factor. With the trace-free conformal background 
extrinsic curvature defined by Aj = ip 6 (Kj — \5]K) and 



with K — 0, we can rewrite equation (Al) as 



n 



(fe) 



1 

8^ 



^-iC(k)d 2 Sj. 



(A2) 



This form of the equation has the advantage that it does 
not involve the conformal factor ip, so we can compute the 
momentum without having a solution of the Hamiltonian 
constraint. Also, I emphasize that we have not used the 
fact that ij) — > 1 at spatial infinity in deriving equation 
( A2 ) , so that if £^ represents an exact symmetry of the 



physical metric, then equation ( A2) can be evaluated (for 
that Killing vector) by integrating over any two-surface 
containing the support of the gravitational field. 

The concept of the momenta (linear and angular) of 
an individual hole in the presence of other holes cannot 
be rigorously defined in general relativity. However, a 
reasonable quasi-local definition for th e momenta of an 
individual black hole is equation (A2) integrated over 
a two-surface exterior to that hole. For evaluating the 
components of the hole's linear momentum, we use the 
three translational Killing vectors of flat Euclidean 3- 
space, and for the angular momentum, we use the three 
rotational Killing vectors with the origin of rotation cho- 
sen to be the center of the hole. 

Following the conformal-imaging approach, the back- 
ground extrinsic curvature of a spatial slice is constructed 
as the linear sum of "single-hole" extrinsic curvature so- 
lutions plus image terms that maintain an isometry con- 
dition (cf Cook 0] ) • The background extrinsic curvature 
for a single black hole (including self-image terms) is pa- 
rameterized directly in terms of the physical momenta 
measured at infinity. Evaluation of the black hole's mo- 
menta via the quasi-local definition is independent of the 
radius of the surface on which the integral is evaluated 
and always yields the correct physical result. Now con- 
sider evaluating the quasi-local momentum integral over 
a two-surface that does not contain the black hole. Using 
Gauss' law, we can rewrite the integral as 



<PV = 0, 



(A3) 



since £^ is a Killing vector of 7^ and A\ satisfies the 
vacuum momentum constraint equation in this volume. 
Constructing a multi-hole extrinsic curvature from single- 
hole extrinsic curvature solutions, including self-image 
terms but not including general image terms, we see that 



equation (A3) implies that the contributions to the ex- 
trinsic curvature from additional holes do not effect the 
quasi-local momenta of a given hole. 

What remains is to examine the contribution of gen- 
eral image terms to the quasi-local momenta of a hole. It 
seems reasonable that these terms should have no contri- 
bution, however, I have so far been unable to prove this 
analytically. Fortunately, it is straightforward to show 
numerically that general image terms make no contribu- 
tion to the quasi-local momenta of either hole in a general 
binary configuration. More specifically, we can construct 
a general solution for the background extrinsic curvature 
including any number of image terms. Computing the 
quasi-local momenta for either of the holes, we can use 
Richardson extrapolation to show that, up to the numer- 
ical precision of the computer, the results are identical 
to those obtained in the absence of any general image 
terms. This is independent of the sizes and separations 
of the holes and of the number of image terms included 
in the extrinsic curvature implying that each image term 
independently contributes nothing to any of the surface 
integrals. I suspect that this result holds for any number 
of holes, however, this has not been verified. 

We see then that within the limitations of defining the 
momenta of an individual hole, the momenta used to 
parameterize a single hole are the momenta of individual 
holes within a multi-hole configuration. That is, linear or 
angular momenta on one hole do not induce any amount 
of linear or angular momentum on any other hole in the 
system. In particular, the orbital angular momentum of 
a system of holes is well defined and does not contain 
an induced spin on any of the holes due to the linear 
momenta of the holes. 
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